International  Journal  of  1 


:  and  Mass  Transfer  80  (2015)  227-235 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 


International  Journal  of  Heat  and  Mass  Transfer 


epage:  www.el 


Characteristics  analysis  and  parametric  study  of  a  thermoelectric 
generator  by  considering  variable  material  properties  and  heat  losses 


CrossMark 


Jing-Hui  Menga,  Xin-Xin  Zhang a,  Xiao-Dong  Wangb,c’* 

a  School  of  Mechanical  Engineering,  University  of  Science  and  Technology  Beijing,  Beijing  100083,  China 

b State  Key  Laboratory  of  Alternate  Electrical  Power  System  with  Renewable  Energy  Sources,  North  China  Electric  Power  University,  Beijing  102206,  China 
c  Beijing  Key  Laboratory  of  Multiphase  Flow  and  Heat  Transfer  for  Low  Grade  Energy,  North  China  Electric  Power  University,  Beijing  102206,  China 


ARTICLE 


INFO 


A  B  S  T  R 


C  T 


Article  history: 

Received  30  November  2013 
Received  in  revised  form  23  July  2014 
Accepted  11  September  2014 
Available  online  29  September  2014 


Keywords 

Thermoelectric  generator 
Variable  properties 
Heat  losses 
Optimization 
Output  power 
Conversion  efficiency 


The  output  power  and  conversion  efficiency  of  the  thermoelectric  generator  (TEG)  are  closely  related  to 
not  only  the  materials  properties  but  also  the  geometric  structure.  This  paper  developed  a  multi-physics, 
steady-state,  and  three-dimensional  numerical  TEG  model  to  investigate  the  TEG  performance,  and  then 
the  model  is  compared  with  the  classical  thermal  resistance  model.  Bismuth-telluride  are  used  as  p-  and 
n-type  materials.  The  comparison  reveals  that  the  assumption  of  constant  material  properties  leads  to 
underestimated  inner  electrical  resistance,  and  overestimated  thermal  conductance  and  Seebeck  coeffi¬ 
cient,  so  that  the  thermal  resistance  model  predicts  unrealistically  high  performance  than  the  present 
model.  The  results  also  indicate  that  when  heat  losses  exist  between  the  TEG  and  the  ambient,  although 
the  output  power  is  slightly  elevated,  the  conversion  efficiency  is  significantly  reduced,  hence,  improve¬ 
ment  of  the  heat  insulation  effect  is  critically  important  for  high-temperature  TEGs.  Furthermore,  the  TEG 
geometry  also  affects  its  performance  significantly:  usage  of  thin  ceramic  plates  increases  the  junction 
temperature  difference,  and  hence  enhances  the  TEG  performance;  there  are  two  optimal  leg  lengths 
which  correspond  to  the  maximum  output  power  and  the  maximum  conversion  efficiency,  respectively: 
when  heat  losses  are  not  ignorable,  a  large  semiconductor  cross-sectional  area  remarkably  reduces  the 
ratio  of  the  heat  liberated  to  the  ambient  to  the  heat  absorbed  from  the  high-temperature  heat  source, 
and  hence  improves  the  conversion  efficiency. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Energy  conservation  and  emissions  reduction  issue  has  become 
a  global  consensus  due  to  the  serious  greenhouse  effect  and 
climate  change.  Hence,  renewable  energy  sources,  such  as  wind, 
solar,  hydro  and  geothermal  energy,  have  made  a  considerable 
development  and  progress  in  the  last  decades.  As  one  of  such  alter¬ 
natives,  thermoelectric  generator  (TEG)  directly  converts  heat  into 
electricity  by  Seebeck  effect  of  thermoelectric  materials.  The  TEG  is 
more  robust,  flexible  and  reliable  when  compared  with  convec- 
tional  power  generators,  because  it  requires  neither  real  working 
fluids  nor  moving  parts,  so  that  it  is  attracting  more  and  more 
attention  [1-12]. 

The  output  power  and  conversion  efficiency  as  indicators  for 
evaluating  the  TEG  performance  are  closely  related  to  the 
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figure-of-merit  of  thermoelectric  materials,  ZT  =  a20T/l,  where  a 
is  the  Seebeck  coefficient,  a  is  the  electric  conductivity,  X  is  the 
thermal  conductivity,  and  T  is  the  absolute  temperature  at  which 
the  properties  are  measured  [13-15],  Thermoelectric  materials 
with  larger  ZT  can  improve  the  TEG  performance.  Previous  research 
has  demonstrated  that  a,  a,  and  X  for  most  thermoelectric  materi¬ 
als  are  all  strongly  temperature-dependent.  It  should  be  noted  that 
the  TEG  generally  works  at  a  higher  temperature  difference  than 
the  thermoelectric  cooler  (TEC),  such  as  the  recycling  of  waste  heat 
for  iron  and  steel  industry.  Thus,  it  is  expected  that  the  tempera¬ 
ture-dependent  materials  properties  have  much  significant  effect 
on  the  TEG  performance.  In  the  previous  research,  the  thermal 
resistance  model  was  widely  used  to  investigate  the  TEG  perfor¬ 
mance,  and  the  model  is  still  most  frequently  adopted  so  far.  This 
kind  of  model  can  derive  analytical  expressions  of  the  output 
power  and  conversion  efficiency,  so  that  it  can  be  applied  for  pre¬ 
liminary  screening  of  TEG  design  and  performance  estimation. 
However,  the  model  accounts  for  energy  balances  only  at  the  hot 
end  and  the  cold  end  of  the  TEG.  As  a  result,  the  constant 
[1-4,6,8,10-12],  or  temperature-averaged  [5]  materials  properties 
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must  be  assumed.  Recently.  Chen  et  al.  [9]  proposed  and  imple¬ 
mented  a  3D  numerical  model  for  TEG  in  FLUENT  UDS  (User 
Defined  Scalar)  environment.  Using  the  same  model,  Reddy  et  al. 
[16]  investigated  the  performance  of  a  novel  composite  thermo¬ 
electric  device  where  a  major  part  of  the  semiconductor  is  replaced 
with  a  conductor  such  as  copper.  In  their  studies  [9,16],  the  great 
influence  of  materials  properties  are  emphasized  on  the  TEG  per¬ 
formance.  We  also  developed  general,  steady-state  or  transient, 
3D  numerical  models  of  thermoelectric  devices  [17-19],  in  which 
the  heat  conduction  equation  and  the  electric  potential  equation 
were  solved  coupled  and  all  thermoelectric  effects  were  taken  into 
account,  including  Seebeck  effect,  Peltier  effect,  Thomson  effect, 
Joule  heating,  Fourier’s  heat  conduction,  and  heat  loss  to  the  ambi¬ 
ent.  The  materials  properties  were  proven  to  have  significant  effect 
not  only  on  the  steady-state  performance  but  also  on  dynamic 
characteristics  for  TEC.  However,  there  is  no  detailed  explanations 
in  the  literature  how  materials  properties  affect  the  TEG 
performance. 

In  addition  to  the  materials  properties  (a,  a,  X,  or  ZT),  the 
geometric  structure  of  the  TEG  also  significantly  affects  its  perfor¬ 
mance.  By  optimizing  the  ceramic  plate  thickness,  the  leg  length 
and  the  cross-sectional  area  of  semiconductors,  a  high-performance 
TEG  can  be  designed.  Some  research  has  already  studied  the  geo¬ 
metric-dependent  TEG  performance  using  the  thermal  resistance 
model  [20-22],  or  one-dimensional  model  [23,24],  However,  the 
three-dimensional  temperature  effect  cannot  be  ignored  when 
there  exist  significant  heat  losses  from  the  TEG  to  the  ambient. 
Recently,  Jang  et  al.  [25]  investigated  the  optimal  structure  of 
high-performance  micro-TEG  at  room  temperature  using  a  three- 
dimensional  model.  However,  constant  materials  properties  are 
adopted  in  their  work,  because  the  TEG  operates  with  only  a  15  K 
temperature  difference. 

This  work  develops  a  multi-physics,  steady-state,  and  three- 
dimensional  numerical  TEG  model  to  investigate  the  optimal 
TEG  geometric  structure.  The  proposed  model  is  compared  with 
the  classical  thermal  resistance  model  to  demonstrate  the  effects 
of  temperature-dependent  materials  properties  and  heat  losses 
on  the  accurate  prediction  of  the  TEG  performance.  Especially,  this 
paper  will  explain  from  the  heat  transfer  mechanism  why  the 
thermal  resistance  model  predicts  unrealistically  high  perfor¬ 
mance  of  TEG.  The  article  structure  arrangement  is  organized  as 
follows.  Sections  2  and  3  describe  and  validate  the  new  three- 
dimensional  TEG  model.  Section  4  compares  the  new  model  with 
the  conventional  thermal  resistance  model  to  highlight  the  effect 
of  temperature-dependent  materials  properties.  Section  5  reveals 
the  reduction  in  the  TEG  performance  due  to  its  heat  losses  to 
the  ambient.  Finally,  the  optimal  structure  parameters  of  the 
TEG  are  discussed  in  Section  6,  including  the  ceramic  plate 
thickness  (Section  6.1),  the  leg  length  (Section  6.2),  and  the 
cross-sectional  area  (Section  6.3)  of  semiconductors. 

2.  Model 

2.1.  Numerical  model 

As  the  basic  element  of  a  practical  TEG,  the  p-n  junction  is 
formed  by  connecting  two  different  types  of  thermoelectric 
materials  (the  hole-rich  p-type  semiconductor  and  the  electron- 
rich  n-type  semiconductor)  together.  When  a  temperature 
difference  is  supplied  to  the  two  sides  of  the  junction,  the  thermal 
excitation  causes  migration  of  carriers  from  the  hot  side  to  the  cold 
side,  so  that  a  Thomson  electromotive  force  will  be  created.  In 
addition,  there  exists  a  concentration  difference  of  carriers 
between  p-  and  n-type  materials,  which  activates  the  electrons  dif¬ 
fusion  from  the  n-type  material  to  the  p-type  material  and  the  back 
diffusion  of  holes.  As  a  result,  a  Peltier  electromotive  force  will  be 


constructed.  The  Seebeck  electromotive  force  is  the  sum  of  the 
Thomson  and  Peltier  electromotive  forces,  which  is  a  maximum 
electric  potential  produced  by  the  junction  and  is  also  referred  as 
to  the  open-circuit  voltage.  However,  when  a  load  is  connected 
to  the  junction,  the  junction  cannot  produce  the  electric  potential 
difference  as  large  as  the  open-circuit  voltage  because  of  the  elec¬ 
tric  potential  drop  caused  by  inner  electrical  resistance  of  the  junc¬ 
tion.  To  obtain  a  higher  output  voltage,  many  p-n  junctions  need  to 
be  connected  in  series  to  form  a  TEG. 

Fig.  1  shows  the  schematic  of  a  TEG  unit,  where  a  p-n  junction 
(n-doped  and  p-doped  bismuth-telluride  legs)  is  connected  electri¬ 
cally  in  series  by  copper  connectors,  and  is  sandwiched  between 
two  same  ceramic  plates  (silicon).  As  electrical  insulation  layer, 
the  ceramic  plate  plays  a  role  in  thermal  conduction  from  outer 
heat  sources  and  cooling  parts  to  p-n  junction. 

The  three-dimensional,  steady-state  TEG  model  developed  here 
solves  heat  and  electricity  conductions  simultaneously.  The 
governing  equations  and  boundary  conditions  are  as  follows. 

2.1.1.  Heat  conduction  equations 

V  •  (2,-VT)  —  VT  =  0  (1) 

where  X  is  the  thermal  conductivity,  a  is  the  electric  conductivity, 
and  ft  is  the  Thomson  coefficient.  The  subscript  i  denotes  conn  for 
the  connectors,  p  for  the  p-type  semiconductor,  n  for  the  n-type 
semiconductor,  cer  for  the  ceramic  plates,  respectively./ is  the  local 
current  density  vector  and  equals  zero  for  the  ceramic  plates.  The 
first  term  on  the  left  side  in  Eq.  ( 1 )  denotes  the  Fourier  heat  conduc¬ 
tion,  and  the  second  and  third  terms  denote  the  Joule  heating  and 
Thomson  effect,  respectively. 

2.1.2.  Electric  potential  equation 

V  •  (a(V<j>  -  aVT))  =  0  (2) 

where  a  is  the  Seebeck  coefficient  of  semiconductors,  </>  is  the 
electric  potential,  aVT  is  Seebeck  electromotive  force  coming  from 
the  Seebeck  effect. 

The  Thomson  coefficient  is  proportional  to  the  first  derivative  of 
the  Seebeck  coefficient  versus  temperature,  or: 


Once  the  electric  potential  is  obtained,  the  current  density  vector 
can  be  calculated  by  the  following  equation: 

J  =  oE  =  a(- V(j>  +  aVT)  (4) 

2.1.3.  Boundary  conditions 

The  boundary  conditions  are  described  in  Fig.  1.  The  tempera¬ 
tures  at  the  top  of  the  upper  ceramic  plate  and  the  bottom  of  the 
lower  ceramic  plate  are  450  K  and  300  K,  respectively.  At  the 
interfaces  between  different  materials,  temperature  and  heat  flux 
are  all  assumed  continuous.  The  side  surfaces  of  n-  and  p-type 
semiconductors  are  convective  boundary  conditions,  or: 

-ig-V-U  (5) 

where  h  is  the  convective  heat  transfer  coefficient,  T* ,  is  the 
ambient  temperature  assumed  to  be  300  K.  The  boundary  condi¬ 
tions  for  electric  conduction  include:  (1)  at  the  side  surface  of  the 
connector  under  the  n-type  semiconductor:  J  =  I/Ac  =  constant, 
where  I  is  the  load  current  and  Ac  is  the  side  surface  area  of  the 
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Hot  end,  7  H=450K 


connector;  (2)  at  the  side  surface  of  the  connector  under  the  p-type  correlations  in  Table  1 ;  (8)  and  then  return  to  step  (4)  until  the 
semiconductor:  </>  =  0.  convergent  iterative  solutions  are  achieved. 


2.1.4.  The  TEG  performance 

Generally,  the  performance  of  the  TEG  is  evaluated  by  two 
parameters,  the  output  power  P  and  conversion  efficiency  t],  which 
are  defined  as  follows,  respectively: 


P  =  N  (6) 


where  V  is  the  output  voltage,  Qh  is  the  heat  supplied  to  the  hot  end 
of  the  TEG. 


2.2.  Solution  procedure 

The  solution  procedure  for  the  present  coupled  model  is 
described  briefly  in  the  following:  (1)  make  the  initial  guess  for 
temperature  and  electric  potential;  (2)  create  geometry  and  grids 
of  the  TEG  model;  (3)  specify  all  boundary  conditions;  (4)  solve 
the  electric  potential  by  Eq.  (2);  (5)  calculate  the  local  current 
density  vector  by  Eq.  (4);  (6)  solve  the  temperature  by  Eq.  (1); 
(7)  calculate  the  Seebeck  coefficient,  Thomson  coefficient,  electric 
conductivity,  and  thermal  conductivity  of  all  materials  by  the 


2.3.  Thermal  resistance  model 


One  of  the  objectives  of  this  work  is  to  compare  the  new  three- 
dimensional  TEG  model  with  the  conventional  thermal  resistance 
model,  and  hence  the  thermal  resistance  model  is  also  introduced 
briefly  here.  Restated  that  the  thermal  resistance  model  does  not 
solve  the  distributions  of  temperature  and  electric  potential  within 
the  TEG,  and  it  considers  heat  balances  only  at  the  hot  and  cold 
ends  of  the  TEG.  Consequently,  the  total  electrical  resistance  R,  See¬ 
beck  coefficient  A  and  thermal  conductance  K  for  the  p-n  junction 
must  be  defined,  or: 


tfp 

rSp 


Hn 

<T„Sn 


A  =  OCp  -  On 


Vp 
-  H 


2nSn 

Hn 


(8) 

(9) 


(10) 


where  Hp  =  Hn  =  H2,  Sp  =  Sn  =  L2x  l2. 

The  heat  balance  equations  for  the  hot  and  cold  ends  can  be 
respectively  expressed  as: 


Table  1 

Material  properties  used  in  the  present  simulations. 


Parameters 

li(wr' m-1) 

p  (£1  m) 
a  (V IT1) 


Semiconductor 


Constant  property 


1.54 
Pp  =  Pn 
1.03  x  10-5 

2.0  x  105 
[26] 


Variable  property 


0.000029T2  -  0.019593T  +  4.809677 
Pp  =  P » 

10-6(0.043542T  -  2.754139) 

10-6(— 0.002025T2  +  1.423448T  -  44.953611) 
[26] 


Connector  (Cu) 


350 

1.695  x  10  9 
6.5  x  10-6 


[25] 


Ceramic  plate  (Si) 


130 


[25] 


230 
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Qh  =  (Th  -  rHJ)D„  (11) 

Qh  =  AITh j  -^/2R  +  K(THj  -  rCJ)  (12) 

Qc  =  A/TCJ  +1/2R  +  /f(r„j  -  Tcj)  (13) 

Qc  =  (Tcj-Tc)Dc  (14) 

P=  Qh  -  Qc  =  A/(THj  -  rcj)  -  I2R  (15) 


where  Qh.  I  Th.  Tc  and  P  have  the  same  meaning  as  in  the  numerical 
model,  Qc  is  the  heat  flow  rate  liberated  from  the  TEG  unit  to  the 
heat  sink,  THj  and  TCj  are  the  temperatures  of  the  hot  and  cold  junc¬ 
tions,  Dh  and  Dc  are  the  thermal  conductances  of  the  upper  and 
lower  ceramic  plates.  The  heat  flows  AITh  j  and  AITCj  caused  by 
the  Peltier  effect  are  released  from  the  hot  junction  and  absorbed 
at  the  cold  junction,  respectively.  Joule  heating  I2R  due  to  the  flow 
of  electrical  current  through  the  p-n  junction  is  generated  inside 
the  p-n  junction.  It  should  be  noted  that  the  Joule  heating  in  the 
thermal  resistance  model  is  generally  assumed  to  be  equally  dis¬ 
tributed  to  the  hot  and  cold  junctions.  The  heat  flow  (TH  -  THj)DH 
through  the  ceramic  plate  at  the  hot  end,  (TCj  -  TC)DH  through  the 
ceramic  plate  at  the  cold  end,  and  K(TH j  —  TCj)  through  the  p-n 
junction  are  caused  by  the  Fourier  heat  conductions. 

3.  Model  validation 

The  same  bismuth-telluride  materials  with  temperature- 
dependent  properties  are  adopted  as  the  p-type  and  n-type  semi¬ 
conductors.  The  silicon  and  copper  are  respectively  selected  as 
the  ceramic  plate  and  connector,  their  properties  are  assumed  to 
be  constant  due  to  ignorable  thermoelectric  effect.  The  properties 
for  all  the  materials  can  be  found  in  Table  1. 

The  grid  sensitivity  analysis  is  performed  for  the  present  TEG 
model  to  check  its  grid  independence,  where  four  grid  systems 
are  designated  as  grid  (i),  grid  (ii),  grid  (iii)  and  grid  (iv), 
respectively.  The  geometric  parameters  of  the  TEG  unit  and  corre¬ 
sponding  grid  number  are  listed  in  Table  2.  The  temperature  differ¬ 
ence  between  the  hot  and  cold  ends  is  assumed  to  be 
AT  =  Th  -  Tc  =  150  K  with  the  convective  heat  transfer  coefficient 
h  =  0  W  m  2  K  '.  Fig.  2  shows  P  and  r}  of  the  TEG  unit  as  functions 
of  I  for  the  four  grid  schemes.  At  /  <  0.1  A,  there  are  almost  no  dif¬ 
ferences  in  P  and  t]  predicted  by  the  four  grid  schemes,  while 
I  >  0.1  A,  the  differences  for  P  and  t]  occur,  and  become  more  obvi¬ 
ous  as  the  current  growing.  The  maximum  deviation  between  grids 

(ii)  and  (iii)  is  1.07%,  it  reduced  to  0.64%  for  grids  (iii)  and  (iv).  Grid 

(iii)  was  chosen  for  the  simulations  as  a  tradeoff  between  accuracy 
and  execution  time. 

The  present  TEG  model  is  compared  with  one-dimensional 
steady-state  analytical  solution  of  Angrist  et  al.  [27]  and  numerical 
model  of  Chen  et  al.  [9],  The  geometry  and  material  properties  used 
here  are  the  same  with  Angrist  et  al.’s  solution  and  Chen  et  al.’s 
simulation.  The  TEG  device  includes  127  units,  each  unit  has  the 
geometry  of  Hj  =  0.2  mm,  H2  =  1.6mm,  i,  =  0.5  mm,  and 
L2  =  1.4mm,  as  well  as  constant  material  properties  of  ap  = 
-an  =  226.8  x  10  6VK  \  kp  =  kn  =  1.52  Wm1  K1  and  <rp  = 
<7n  =  1.447  x  10  5 am.  The  TEG  operates  at  TH  =  423K  and 


Fig.  2.  Grid  independence  examination. 


Comparisons  of  numerical  and  analytical  results  (TH  =  423  K  and  Tc  =  303  K). 
Quantity  Present  results  Analytical  results  [27]  Chen  et  al.’s  results  [9] 
Qh  (W)  81.2  81.3  81.8 

P  (W)  3.88  3.98  3.62 

1(A)  1.08  1.08  1.032 

n  4.77  4.89  4.43 


Tc  =  303  K,  and  the  heat  losses  to  the  ambient  are  ignored.  With 
these  conditions,  Qh,  P,  I,  and  »;  are  predicted  by  the  three  models. 

Table  3  shows  that  the  present  prediction  has  better  agreement 

with  the  analytical  solution  than  Chen  et  al.’s.  It  should  be  pointed 
that  the  analytical  solution  [27]  is  strictly  true  only  when  the  TEG 
has  one-dimensional  temperature  distribution,  uniform  current 
density,  constant  material  properties,  and  no  heat  losses  to  the 
ambient.  However,  the  present  model  is  more  general  with 
three-dimensional  distributions  of  temperature  and  current  den¬ 
sity,  which  can  be  used  to  study  the  TEG  performance  at  more 
complicated  operating  conditions  such  as  high  temperature  differ¬ 
ence  and/or  large  heat  losses  to  the  ambient. 


4.  Comparison  between  three-dimensional  model  and  thermal 

resistance  model 

As  stated  above,  the  thermal  resistance  model  is  the  most  fre¬ 
quently  adopted  in  literature,  and  the  model  can  only  use  constant 
materials  properties  because  the  inner  temperature  distribution  of 
the  TEG  is  not  solved.  To  improve  the  accuracy,  the  reference 
temperature  for  materials  properties  is  generally  taken  as 
Tref  =  (Th  +  Tc)/2.  However,  the  present  model  uses  temperature- 
dependent  properties,  thus,  comparison  between  the  two  models 
can  highlight  the  effect  of  temperature-dependent  materials 
properties.  The  comparison  is  performed  at  AT  =  30  K,  60  K,  90  K, 
120  K  or  150  K  with  a  fixed  Tc  =  300  K. 


The  geometric  parameters  and  test  grid  numbers  for  TEG. 


Geometry  Value  t,  I2  t3  H,  ff2  ff3 


Grid  (i)  Grid  number  2  10  4  4  2  20 

Grid  (ii)  Grid  number  3  15  6  6  3  30 

Grid  (iii)  Grid  number  4  20  8  8  4  40 

Grid  (iv)  Grid  number  5  25  10  10  5  50 
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Fig.  3  shows  the  I-P  curves  for  both  the  models  at  different  AT. 
The  I-P  curves  are  all  typical  parabola  for  the  two  models,  P  firstly 
increases  and  then  drops  as  the  current  increases,  hence  there 
exists  a  maximum  output  power,  Pmax.  Pmax  is  elevated  when  AT 
is  increased.  In  addition,  both  the  models  predict  that  as  AT  is 
increased,  the  effective  working  currents,  defined  as  the  current 
when  P  0,  are  all  increased,  for  example,  at  AT  =  30  K,  the  effec¬ 
tive  working  current  interval  for  the  present  model  is  0  <  I  <  0.13  A, 
and  increases  to  0  <  I  <  0.23  A  when  AT  =  60  K. 

Fig.  3  also  indicates  that  although  the  variation  trend  of  P  as  the 
function  of  /  and  AT  predicted  by  the  thermal  resistance  model  is 
similar  with  the  present  model,  the  quantitative  differences 
between  the  two  models  are  significant,  especially  for  high  I  and 
large  AT.  The  differences  include  that  the  thermal  resistance  model 
overestimates  P  and  the  maximum  effective  working  current.  The 
overestimation  can  be  explained  as  follows.  Because  the  present 
model  has  solved  distributions  of  the  temperature  and  electrical 
potential,  the  open-circuit  voltage  (V0)  and  the  output  voltage 
can  be  directly  obtained.  Thus,  the  effective  electrical  resistance 
can  be  calculated  as  R  =  (Va  -  V)/I,  which  is  shown  in  Fig.  4.  It  can 
be  seen  that  R  in  the  thermal  resistance  model  remains  constant 
according  to  Eq.  (8),  however,  R  in  the  present  model  is  higher  than 
that  in  the  thermal  resistance  model,  and  their  difference  becomes 
larger  as  I  increases.  Consequently,  the  overestimation  of  P  in  the 
thermal  resistance  model  can  be  attributed  to  the  underestimation 


Current  (A) 


Fig.  3.  The  output  power  of  the  TEG  unit  for  constant  and  variable  properties  at 
various  temperature  differences. 


of  R  according  to  Eq.  (15).  In  addition,  with  P  =  0,  Eq.  (15)  can  be 
rearranged  as  follows: 

(16) 

where  /max  is  the  maximum  working  current.  Based  on  the  same 
reason,  fmax  in  the  thermal  resistance  model  is  also  overestimated. 

Fig.  5  shows  the  heat  absorbed  at  the  hot  end  of  the  TEG  unit  for 
the  two  models  at  various  AT.  Compared  with  the  present  model, 
the  thermal  resistance  model  predicts  a  larger  Qh,  especially  for 
large  /.  According  to  Eq.  (12),  Qh  can  be  split  into  three  parts:  the 
Peltier  heat,  Joule  heat  and  Fourier  heat  conduction,  their  contribu¬ 
tions  for  AT  =  1 50  K  are  presented  in  Fig.  6.  Due  to  the  quadratic 
relationship  between  the  Joule  heat  and  I,  the  Joule  heat  calculated 
by  the  two  models  has  only  a  small  difference.  As  a  plain  heat 
source,  Peltier  heat  is  generated  only  on  the  interface  between 
two  semiconductors  or  two  conductors.  The  reference  temperature 
for  a (T)  in  the  Peltier  heat  (=a (T)IT)  should  be  taken  as  the  interface 
temperature.  Thus,  we  have  a  =  a(THj)  for  the  hot  junction.  How¬ 
ever,  when  the  approach  with  temperature-averaged  properties 
are  adopted,  almost  all  of  thermal  resistance  models  assume 
a  =  a(Tref)  [22,28-31],  For  bismuth-telluride  materials,  a  has  a 
negative  temperature  effect,  a(Tref)  is  larger  than  a(THj)  due  to 


Current  (A) 


Fig.  5.  The  heat  absorbed  at  the  hot  end  of  the  TEG  unit  for  constant  and  variable 
properties  at  various  temperature  differences. 
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Fig.  4.  Variations  of  the  total  electrical  resistance  and  the  total  thermal  conduc¬ 
tance  with  /  at  AT  -  150  K. 
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Fig.  6.  The  produced  heat  within  the  TEG  unit  for  AT  =  150  K. 
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properties  at  various  temperature  differences. 


7ref<  ThJ.  Therefore,  the  thermal  resistance  model  overestimates 
the  Peltier  heat  on  the  hot  junction,  as  shown  in  Fig.  6.  According 
to  Eq.  (10),  the  total  thermal  conductance  of  the  p-n  junction  is 
7.70  x  10~4  W  K  1  for  the  thermal  resistance  model.  However,  this 
thermal  conductance  is  actually  related  to  AT  and  the  functional 
relation  between  X  and  T  of  p-  and  n-type  semiconductors,  which 
can  be  calculated  by  solving  only  the  heat  conduction  equation 
(Eq.  (1))  for  the  same  TEG  geometry  and  the  same  AT,  while 
all  thermoelectric  effects  are  ignored.  With  AT=150K,  the 
thermal  conductance  for  the  p-n  junction  is  found  to  be  only 
4.08  x  10~4  W  K  1  (Fig.  4).  Thus,  the  thermal  resistance  model 
overestimates  the  contribution  of  Fourier  heat  conduction.  As  a 
result,  the  thermal  resistance  model  predicts  a  larger  Qh. 

Fig.  7  shows  r\  of  the  TEG  unit  as  a  function  of  I  at  various  AT  for 
the  two  models.  Comparison  of  Figs.  3  and  5  indicates  that 
although  P  and  Qh  are  all  overestimated  by  the  thermal  resistance 
model,  the  overestimation  of  P  is  more  significant.  Hence,  r\  is  also 
overestimated  significantly  by  the  thermal  resistance  model.  For 
example,  the  maximum  conversion  efficiencies,  rjmax,  are  3.85% 
and  5.38%  for  the  two  models  at  AT  =150  K,  respectively,  and  the 
thermal  resistance  model  overestimates  r/max  by  39.74%. 

5.  Effect  of  heat  losses 

In  practical  applications,  the  TEG  has  the  promising  potential 
for  the  recycling  of  industrial  waste  heat  such  as  in  iron  and  steel 
industry,  where  the  waste  heat  temperature  is  generally  higher 
than  450  K  [32],  so  that  the  performance  of  the  TEG  is  inevitably 
affected  by  heat  losses  between  the  TEG  and  the  ambient.  It  should 
be  noted  that  the  thermal  resistance  model  cannot  take  heat  losses 
into  account,  hence,  the  effects  of  heat  losses  are  investigated  only 
by  the  present  model  with  temperature-dependent  properties  at 
AT  =150  K.  The  value  of  total  heat  transfer  coefficient  (h)  is 
assumed  to  be  0-200  Wm~2  K~'[  17],  which  combines  contribu¬ 
tions  of  the  convection  and  radiation. 

Fig.  8  shows  variations  of  P  and  rj  for  various  h  at  AT  =150  K.  At 
/  <  0.1 2  A,  h  has  a  ignorable  effect  on  P,  while  at  /  >  0.1 2  A,  its  effect 
is  elevated  gradually  as  I  increases.  The  occurrence  of  heat  losses 
reduces  the  inner  temperature  of  the  TEG,  and  the  reduction  is 
found  to  be  more  significant  for  high  h  (Fig.  9).  For  bismuth-tellu- 
ride  materials,  the  electrical  resistivity  has  a  positive  temperature 
effect,  as  a  result  of  the  reduced  inner  temperature,  the  effective 
electrical  resistance  and  inner  electric  potential  drop  of  the  p-n 
junction  are  both  reduced  when  h  is  increased.  On  the  other  hand, 
the  open-circuit  voltage  remains  almost  unchanged  for  various  h. 


Current  (A) 


Fig.  8.  Performance  of  the  TEG  unit  with  various  h  at  AT  =  150  K. 


Fig.  9.  Temperature  distribution  of  z  direction  along  the  center  of  p-type  semicon¬ 
ductor  with  various  h  at  AT  =  150  K. 

Hence,  the  higher  output  voltage  (V=V0-IR)  and  output  power 
(P  =  IV)  occur  at  a  larger  h.  For  the  same  reason  in  Section  4,  the 
effective  working  current  interval  increases  as  h  is  increased.  At 
li  =  0Wm  2 K  1 ,  the  maximum  effective  working  current  is 
0.45  A,  it  increases  to  0.46  A  when  h  =  200  W  m  2  K  \ 

Although  the  heat  losses  slightly  elevate  the  output  power,  the 
conversion  efficiency  is  significantly  reduced  as  h  increases.  A 
higher  h  means  that  more  heat  will  be  transferred  to  the  ambient 


:  AT  =  150  K. 
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from  the  TEG,  hence,  a  larger  Qh  will  be  needed  to  keep  the  con¬ 
stant  AT  =  150  I<  (Fig.  10).  The  result  indicates  that  the  high  output 
power  comes  at  the  cost  of  conversion  efficiency,  thus,  the  heat 
losses  should  be  reduced  as  small  as  possible  when  a  TEG  operates. 
For  example,  as  shown  in  Fig.  10,  Pmax  is  5.74  x  10  3  W  for 
h  =  0  W  m  2  K  1  and  5.85  x  10  3  W  for  h  =  200  W  m2  K  \  which 
is  increased  by  only  1.92%,  however,  t/max  is  decreased  by  19.43% 
(=(3.86%  -  3.1 1%)/3.86%). 

6.  Parametric  study  of  TEG  geometry 

Effect  of  the  TEG  geometry  on  the  performance  is  investigated 
by  the  present  model  and  the  thermal  resistance  model.  Because 
the  performance  predicted  by  the  two  models  has  almost  no  qual¬ 
itative  differences  but  quantitative  differences,  and  Section  4  has 
already  explained  why  the  thermal  resistance  model  overestimates 
the  performance  of  the  TEG  unit,  thus,  this  section  will  just  demon¬ 
strate  their  difference  but  not  repeat  the  explanation.  The  output 
power  and  conversion  efficiency  for  the  TEG  with  various 
geometries  are  compared  at  three  characteristic  currents  (small, 
intermediate,  and  large)  for  Sections  6.1  and  6.2  or  two  current 
densities  for  Section  6.3.  The  temperature  difference  is  chosen  as 
AT  =  1 50  K  with  h  =  1 00  W  m  2  K  1  for  the  present  model  and 
h  =  0  W  m  2  K  1  for  the  thermal  resistance  model. 

6.1.  Effect  of  ceramic  plate  thickness 

As  the  electrical  insulator,  the  ceramic  plate  transfers  heat  from 
the  outside  heat  source  to  the  p-n  junction  as  well  as  from  the  p-n 


junction  to  the  outside  heat  sink.  Due  to  non-zero  thermal 
resistance  of  the  ceramic  plate,  its  thickness  directly  affects  the 
temperatures  of  hot  and  cold  junctions,  and  hence  has  significant 
impact  on  the  TEG  performance,  as  shown  in  Fig.  11.  At  the  three 
characteristic  currents  of  0.1,  0.2,  and  0.4  A,  P  and  rj  decrease  line¬ 
arly  as  the  ceramic  plate  thickness  increases,  this  is  mainly  because 
the  temperature  difference  (THj  -  Tc j)  between  the  hot  and  cold 
junctions  becomes  smaller  as  the  plate  gets  thicker,  which  reduces 
the  Seebeck  electromotive  force,  and  hence  reduces  the  open- 
circuit  voltage.  Fig.  1 1  confirms  again  the  difference  in  the  perfor¬ 
mance  predicted  by  two  models,  and  the  difference  becomes  more 
significant  at  high  /,  therefore,  the  thermal  resistance  model  is  only 
applicable  for  preliminary  scanning  of  TEG  performance. 

6.2.  Effect  of  semiconductor  leg  length 

Effect  of  semiconductor  leg  length  on  the  performance  of  the 
TEG  unit  is  shown  in  Fig.  12.  As  the  leg  length  increases,  P  and  >; 
first  reach  their  peak  values  (Pmax  and  rjmBX)  at  a  certain  length, 
then  monotonically  drop.  It  is  worth  noting  that  for  the  three  char¬ 
acteristic  currents,  the  optimal  leg  length  corresponding  to  Pmax  is 
almost  the  same,  however,  the  leg  length  corresponding  to  r/max  is 
different  from  each  other,  which  moves  toward  longer  length 
direction  at  low  current.  Thus,  to  design  high-performance  TEGs, 
one  has  to  make  a  choice  between  the  two  optimal  leg  lengths. 

Fig.  12(a)  also  shows  that  Pmax  is  elevated  as  I  increases,  but  the 
maximum  leg  length  is  reduced.  For  example,  the  maximum  length 
predicted  by  the  present  model  is  larger  than  5.00  mm  at  I  =  0.1  A, 
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however,  the  maximum  length  is  reduced  to  2.40  mm  at  1  =  0.2  A 
and  1.15  mm  at  I  =  0.4  A.  This  can  be  explained  using  the  following 
relation: 

Rm3X  =A(Tu>  ~  Tcj>  (17) 

With  the  fixed  ceramic  plate  thickness  and  semiconductor  materi¬ 
als,  the  open-circuit  voltage,  A(THj  -  TCj),  are  almost  the  same  for 
different  leg  lengths  and  currents.  According  to  Eq.  (17),  the  maxi¬ 
mum  of  the  total  electrical  resistance  of  the  p-n  junction  should 
be  less  for  a  higher  I,  thus,  the  maximum  leg  length  becomes  small. 

6.3.  Effect  of  cross-sectional  area  of  semiconductor 

Effect  of  cross-sectional  area  of  semiconductors  (Sp  =  Sn)  on  the 
TEG  performance  is  investigated  at  TH  =  450  K  and  Tc  =  300  K.  For 
fair  comparison,  the  same  current  density  (/)  for  various  Sp  and 
the  output  power  per  unit  cross-sectional  area  (Parea  =  P/SP)  are 
used  here. 

Fig.  13  confirms  again  that  the  thermal  resistance  model 
overestimates  the  TEG  performance.  Fig.  13(a)  shows  that  Parea 
almost  remains  unchanged  at  various  Sp  for  the  thermal  resistance 
model,  however,  the  present  model  predicts  that  the  TEG  has  bet¬ 
ter  performance  at  a  small  Sp.  According  to  the  thermal  resistance 
model,  Parea  meets: 


1.8 1 — * — ' — * — ' — * — ' — * — ' — * — ' — * — ' — * — ' — * — 

0.0  0.5  1.0  1.5  2.0  2.5  3.0  3.5  4.0 

Cross-sectinal  area  (mm2) 


(18) 


Eq.  (18)  indicates  that  Parea  is  independent  on  Sp.  However,  because 
the  present  model  solves  the  temperature  field,  the  difference  of  the 
temperature  distribution  within  semiconductors  for  various  Sp 
causes  their  inner  electrical  resistance  changes,  and  the  p-n  junc¬ 
tion  with  a  smaller  Sp  has  a  lower  inner  electrical  resistance,  so  that 
a  larger  Parea  occurs  at  the  smaller  Sp. 

Fig.  3b  shows  the  variation  of  tj  with  Sp.  According  to  Eq.  (12), 
the  absorbed  heat  per  unit  cross-sectional  area  (Q^a  =  Qh/Sp)  is  also 
independent  on  Sp,  thus,  r\  =  PlleJ(hrea  remains  almost  unchanged 
due  to  the  same  energy  input  and  output  for  the  thermal  resistance 
model.  However,  for  the  present  model,  rj  increases  first  rapidly 
then  slowly  with  the  increased  Sp.  Restated  that  the  present  model 
takes  heat  losses  (Q]OSS)  into  account,  although  the  TEG  with  a  large 
Sp  has  a  large  total  convective  heat  transfer  area  (Sconv,totai).  the  ratio 
of  convective  heat  transfer  area  to  semiconductor  cross-sectional 
area  (Sconv,totai/Sp)  reduces,  which  implies  that  the  heat  losses  per 
unit  semiconductor  cross-sectional  area  (0iOss,area  =  Qioss/Sp) 
becomes  smaller  for  the  large  Sp.  For  example,  QiOSs,area=  117683 
W  m 1  2  for/  =  8  x  10s  A  m  2 *  and  L2xL2  =  0.25  mm2,  which  reduces 
to  59845  W  m  2  for  /  =  8  x  105  A  m  2  and  I2  x  l2  =  1  mm2.  Conse¬ 
quently,  to  maintain  a  fixed  temperature  difference  1 50  K,  more 
heat  needs  to  be  input  into  the  hot  end  for  the  TEG  with  small  Sp, 
so  that  t]  reduces  as  Sp  decreases.  The  results  above  indicate  that 
for  a  practical  TEG  operation,  when  the  heat  losses  are  not  ignor- 
able,  a  large  Sp  should  be  recommended. 


7.  Conclusions 


This  paper  develops  a  three-dimensional  and  steady-state  TEG 
model,  which  couples  the  temperature  and  electric  potential  equa¬ 
tions  and  accounted  for  all  thermoelectric  effects,  including  the 
Seebeck  effect,  Peltier  effect,  Thomson  effect,  Joule  heating  and 
Fourier’s  heat  conduction.  The  model  is  used  to  investigate  effects 
of  the  temperature-dependent  material  properties,  heat  losses,  and 
TEG  geometry  on  the  TEG  performance.  The  main  conclusions  are 
as  follows: 


(1)  The  present  model  is  compared  with  the  classical  thermal 
resistance  model  in  a  wide  range  of  currents  and  tempera¬ 
ture  differences.  The  comparison  shows  that  the  thermal 
resistance  model  overestimates  the  output  power  and 
conversion  efficiency  of  the  TEG,  it  is  especially  true  at  large 
currents  and  high  temperature  differences.  The  overesti¬ 
mated  performance  can  be  attributed  to  the  fact  that  the 
thermal  resistance  model  does  not  solve  the  temperature 
distributions  within  the  TEG  and  can  only  adopt  the  assump¬ 
tion  of  the  constant  or  temperature-averaged  material  prop¬ 
erties.  For  bismuth-telluride  materials  used  in  this  work,  the 
electrical  resistivity  has  a  positive  temperature  effect  and 
the  Seebeck  coefficient  has  a  negative  temperature  effect. 
Thermal  resistance  model  with  temperature-averaged 
properties  predicts  a  much  lower  effective  inner  electrical 
resistance  of  the  p-n  junction,  and  the  much  higher  Peltier 
heat  on  the  hot  junction  and  Fourier  heat  conduction,  conse¬ 
quently,  overestimates  the  TEG  performance. 

(2)  Although  the  heat  losses  slightly  elevate  the  output  power, 
the  conversion  efficiency  is  significantly  reduced  as  the  total 
heat  transfer  coefficient  increases.  With  higher  total  heat 

transfer  coefficient,  the  hot  end  of  the  TEG  should  absorb 
more  heat  to  maintain  a  fixed  temperature  difference,  which 
leads  to  a  remarkable  reduction  in  conversion  efficiency. 

Therefore,  for  a  real  TEG  operation,  the  heat  losses  should 

be  reduced  as  small  as  possible. 
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(3)  The  TEG  geometric  parameters  have  significant  effect  on  its 
performance.  As  the  ceramic  plate  becomes  thicker,  the  TEG 
performance  deteriorates  due  to  the  reduced  temperature 
difference  between  the  hot  and  cold  junctions.  For  a  specific 
TEG,  there  are  two  optimal  leg  lengths  of  semiconductors 
corresponding  to  the  maximum  output  power  and  maxi¬ 
mum  conversion  efficiency,  respectively.  When  the  heat 
losses  are  not  ignorable,  a  large  cross-sectional  area  should 
be  recommended. 
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